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Abstract We consider the hadroproduction of W gauge 
bosons in their leptonic decay mode. Starting from the 
leading-order expressions, we show that by defining a 
suitable scaling variable the centre-of-mass dependence 
of the cross sections at the LHC energies can be es¬ 
sentially described by a simple power law. The scaling 
exponent is directly linked to the small-cc behaviour of 
parton distribution functions (PDF) which, at the high 
virtualities involved in W production, is largely dictated 
by QCD evolution equations. This entails a particularly 
simple scaling law for the lepton charge asymmetry and 
also predicts that measurements in different collision 
systems (p-p, p-p, p-Pb Pb-Pb) are straightforwardly 
related. The expectations are compared with the exist¬ 
ing data and a very good overall agreement is observed. 
It is shown that the PDF uncertainty in certain cross- 
section ratios between nearby centre-of-mass energies 
can be significantly reduced by taking the ratios at fixed 
value of scaling variable instead of fixed rapidity. 

Keywords Inclusive W production, hadron colliders, 
hard-processes in proton-lead and heavy-ion collisons, 
parton distribution functions 
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1 Introduction 

The production of W gauge bosons in hadronic colli¬ 
sions is a process which is sensitive to practically all as¬ 
pects of Standard Model, from electro-weak couplings 
to QCD dynamics and the non-perturbative parton con¬ 
tent of the hadrons [1] . One of the most precisely mea- 
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sured observables at hadron colliders is the rapidity (y) 
dependence of the lepton charge asymmetry, C^, 


/dy - /dy 
^ ^ da^^ /dy + da^~ /dy’ 


( 1 ) 


where the charged lepton = e, y) originates from 
the leptonic decay of the W boson. This observable is 
a useful probe of proton parton distribution functions 
(PDFs), in particular, to disentangle the flavour depen¬ 
dence [2113] which is not well constrained by the deep 
inelastic scattering^ Today, the charge asymmetry has 
been studied in detail by the CDF |31[S| and DO |71inilSllS] 
experiments in p-p collisions at the Tevatron as well as 
the ATLAS [TOIITT] . CMS jHUS], and LHCb [HlfTK] ex¬ 
periments in p-p collisions at the LHC. While the broad 
features of the experimental data are well captured by 
fixed-order perturbative QCD calculations dllllT], the 
simultaneous reproduction of the DO data in bins of dif¬ 
ferent kinematic cuts is known to pose difficulties [181 

ng. 

The first measurements of W production in p-Pb 
collisions have recently appeared mMm and vari¬ 
ous observables seem to favour the use of EPS09 nu¬ 
clear PDFs (nPDFs) [23] instead of a naive superposi¬ 
tion of free nucleon PDFs (similar conclusion can be 
expected in the case of other sets of nPDFs [211131 
[23]L In addition, these measurements may also help 
to probe, for the first time, the flavour dependence of 
nuclear modifications in quark densities [20] . The pro¬ 
duction of W bosons in heavy-ion collisions is also of 
paramount importance. Measurements by ATLAS m 
and CMS [2H] in Pb-Pb collisions have revealed that the 
production rate approximately scales with the number 


^unless a deuterium target, complicated by possible nuclear 
corrections, is used. 
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of binary nucleon-nucleon collisions. This is in sharp 
contrast to hadronic observables (high-transverse mo¬ 
mentum hadrons [^1501151] and jets 1 which 

are strongly suppressed as compared to p-p collisions. 
Thus, the leptons from W decays are valuable “messen¬ 
gers” from the initial state of heavy-ion collisions and 
could also be used to constrain the nPDFs [MIIMIIST]. 

In this paper, our main focus is on the centre-of- 
mass energy (-^/s) systematics of the production cross 
sections da^ jdy in hadronic collisions and the conse¬ 
quent scaling properties of the lepton charge asymme¬ 
try. First, in Sectionj^ we show how the scaling laws for 
absolute cross sections and charge asymmetries emerge 
from the relatively simple leading-order expressions. In 
Section we then contrast these expectations against 
next-to-leading order (NLO) computations. Section 
presents comparisons with the existing world data from 
LHC and Tevatron experiments as well as demonstrates 
how PDF uncertainties in some ratios of W cross sec¬ 
tions can be suppressed by carefully chosing the rapid¬ 
ity binning. Finally, we summarize our main findings in 
Section [5] 


2 Derivation of the scaling properties 

2.1 Absolute cross sections 


= (pTe“^)/(-\/s — and the momentum ar¬ 

gument xi = {x 2 PTe^)/{x 2 \fs — PTe~y). The Mandel¬ 
stam variables t and u are 


t = —y/SP'YXie U = —^PT'X2e^. 


(3) 


The PDFs are denoted by (withQ^ 

and the sum runs over all flavours i,j such that the elec¬ 
tric charges of the quarks sum up to ±1. Since the 
total width of the W boson is much smaller than its 
mass, Iw ^ AIwj we can make use of a delta-function 
identity e/(x^ + e^) —>■ Tr6{x), as e —?> 0, to perform the 
remaining integral in Eq. Q. We And 


dydpT 


1 


24s ly sin^ 0w 
_ Pt _^ 

Vl - ^Pt/^w Tt 

Y^l — Ap^/M'^ 
— Ap'^/M^ 

l=p yi-4p2,/M^ 


iWw/\v 


(4) 




eq- - 1 - 65 . .±1 




1± 


1± 


)qfHx2) + 






qf{x^)qf^{x^) 


We consider the inclusive production of W bosons in 
high-energy collisions of two hadrons, Hi and H 2 , fol¬ 
lowed by the decay of W to a charged lepton and a 
neutrino. 


Hi -k H 2 ^ W" -k X ^ -k -k X, 


where the momentum arguments of the PDFs are 

^ M^ey 


2pta/s 


•^2 — 


2pT^/s 


l±^l-4p|/M2 


(5) 

( 6 ) 


Hi -k H 2 ^ W+ -k X ^ -k + X. 


At leading order, the production cross section double 
differential in the charged lepton rapidity y and trans¬ 
verse momentum px reads EHUSg, 



_ ( 3 : 13 ^ 2 ) ^ _ 

ys ) {X1X2S - Miff + 

{i+u±i^uf qf^ {xi , Qfqf^ {x 2 , Q^)+ 

{i + uTi± {xi,Q^)qf^ {x2,Q^) , 

where the symbols aem; and Vij refer to the fine- 
structure constant, weak-mixing angle, and elements of 
Cabibbo-Kobayashi-Maskawa matrix, respectively. The 
mass and width of the W boson are denoted by Mw 
and Tw- The lower limit of the X2 integral is given by 



Let us first consider a situation witlj^ y :$> 0, that is, 
x^ < Xi . In terms of a dimensionless variable ^1 (which 
coincides with xf when px —t Myf/2), 


6 = 


Afw 


(7) 


the momentum fractions in Eq. § become 


± Mw - 

.±- "w 


2pTS^l 


1± 




( 8 ) 


(9) 


At sufficiently small x, the sea-quark densities at high 
^ should be reasonably well approximated by 
a power law [40] 


xqq{x,Q^) Ki xqi{x,Q^) Ki Ni X 


( 10 ) 


^For simplicity, y 3> 0 (j/ <SC 0) should be understood as eJ' ^ 1 
(e^ <C 1) In the remainder of the paper. 
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where the exponent a{Q^) > 0 and the normalizations 
Ni should both be almost flavour independent. Such 
a behaviour (though not exactly a power law [IT]) is 
expected by considering the small-a; and large limit 
(the so-called double logarithmic approximation |32|) of 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi parton evo¬ 
lution equations |T5limiT51H5] and it is also consistent 
with the observations in deep inelastic scattering m 
with the dependence of the exponent a{Q'^) be¬ 
ing roughly logarithmic. However, in what follows, the 
“running” of a{Q^) does not directly show up since we 
will always set = M^. For brevity, we will denote 
a = a{Q'^ = M^) from now on. By using the approx¬ 
imation Eq. (10) in Eq. and trading the rapidity 
variable y with ^i, we hnd 


(s,Ci) 


s“x/±(ei,PT,Hi,H2), 2/»0, (11) 


dprd^i 

where /^(^,pt, Hi, H 2 ) is a function that does not de¬ 
pend explicitly on s or y. Since the expression of Eq. 
is peaked at px ~ and the px dependence of the 

probed momentum fractions in Eqs. § is not particu¬ 
larly herce, the x interval spanned by integration over 
px with a typical kinematic cut px ^ 20 GeV remains 


sufficiently narrow such that approximation of Eq. (101 


stays valid. Under these conditions, the scaling law in 
Eq. 0 holds also for px-integrated cross sections. 


do'^ (s,Ci) 

di\ 


xF±(a,Hi,H2), p»0, 


( 12 ) 


where F±(^i,Hi,H 2 ) = / dpx/=^(Ci,PT,Hi)6»(px-PT“). 
In the backward direction with p <C 0, the appropriate 
scaling variable is 


6 = 




(13) 


[T5]l. At small x, the probed x regions will be different 
for two different y/s, see Eq. §, but at large the 
X dependence is almost purely dictated by the DGLAP 
evolution (in our scaling law approximated by a power 
law) and less prone to PDE uncertainties. We will come 
back to this later on in Section lU^ 


2.2 Gharge asymmetries 


Since the ^/s dependence in Eqs. ( |l2[ ) and (14) is com¬ 
pletely in the common prefactor s“, it follows that the 
lepton charge asymmetry Eq. Q should obey a partic¬ 
ularly simple scaling law. 


X-H 1 .H 2 


(s,6) 

where 


E(ei,Hi,H2), 

G(6,Hi,H2), 


2 /> 0 , 

2/<0, 


(15) 


i^(e,Hi,H2) 


f+(e,Hi,H2)-F-(g,Hi,H2) 

F+(e,Hi,H2)+E-(e,Hi,H2)’ 


(16) 


and similarly for G. In other words, at hxed or ^ 2 j 
the charge-asymmetry should be approximately inde¬ 
pendent of the centre-of-mass energy. In fact, here one 
can allow the exponent a to depend also on ^/s and ^ 1^2 
and it is only required that the PDFs are locally well 
approximated by a power law in the relevant region at 
small-cc. 

Another, and also a bit surprising feature of the 
charge asymmetry is that at sufficiently large \y\ it ef¬ 
fectively depends only on the nucleon that is probed at 
large x. This follows from the facts that when \y\ is suffi¬ 
ciently large, either ud or du partonic process eventually 
dominates, and that the light-sea-quark distributions 
are expected to be approximately SU(2) symmetric at 
small X, 


such that 

^^^-^^«s“xG±(6,Hi,H2), 2/«0, (14) 

d^2 

where G^ (C 2 , Hi, H 2 ) is a function that does not depend 
explicitly on s or y. If Hi = H 2 , then F=*=(^i, Hi, H 2 ) = 

g±(6,Hi,H2). 

Here, we emphasize the fact that at fixed ^1 (^ 2 ) 
the X region at which the PDFs of hadron Hi (H 2 ) is 
sampled becomes approximately independent of yfs, see 
Eq. Q. Going to forward (backward) direction pushes 
this region to large x where the parameterization de¬ 
pendence of PDEs may be large. As a consequence, one 
could hope that the PDE uncertainties on cross-sections 
ratios between two different values of ^/s would better 
cancel out if performed at fixed ^12 than at fixed ra¬ 
pidity (as has been done e.g. by LHGb collaboration 


u{x, Q^) « u{x, Q"^) « d{x, Q^) « d{x, Q^), a; <C 1, 

(17) 

and thus symmetric with respect to charge conjugation 
and isospin rotation. For example, one would expect 
that C^’’’(s,^i) « C|’’^(s,^i) at large ^ 1 . In the case 
of nuclei the nPDFs ff^(x,Q‘^) are built from the free 
nucleon PDFs /f™*°''(a;, Q^) and nuclear modification 
factors g [23] 1 

/f"(x, Q^) = Q^) + Q2), 

(18) 

where 

/f ^ ^proton.A^proton^^^ q2)^ 

/“’^(x,Q^) = (20) 
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At small-a; one expects modest shadowing ^ 

1) which, however, should not significantly alter the 
scaling exponent a (particularly at high ^ in¬ 
volved here) and, to a good approximation, the effect of 
shadowing is just a slight overall downward normaliza¬ 
tion in the absolute cross sections which should largely 
disappear in the case of charge asymmetry. In other 
words, we can encapsulate the scaling law for lepton 
charge asymmetry as 

2 /» 0 , 

2 /« 0 , ( 21 ) 

independently of the nature of hadron (nucleon, anti¬ 
nucleon, nucleus) probed at small x. 



3 Scaling vs. NLO calculation 


Most of our plots in the rest of the paper will use the 
scaling variables which are related to rapidity y and 
centre-of-mass energy ^/s via Eq. Q and Eq. (13). To 
ease the interpretation in what follows, this dependence 
is illustrated in Fig. 



Fig. 1 Relation of rapidity y and scaling variables for a 
few values of ^/s. 


According to Eq. (10), the scaling exponent a in 
Eq. (12) should reflect the small-x behaviour of quark 
distributions and it can be straightforwardly extracted 
from cross-sections at two different centre-of-mass ener¬ 
gies. To verify this correspondence and the consistency 
of our derivation, we have computed the full NLO cross- 
sections at ^/s = 7, 8,13TeV for p-p collisions using 
MCFM Monte-Carlo code [35] and CTIONLO PDFs [I^ . 
From these cross-sections, we have evaluated the effec- 



Fig. 2 Scaling exponent extracted from NLO calculations 
(upper panel) and its comparison with CTIONLO PDFs 
(lower panel). 


tive scaling exponent cteff by 


aesiO = log 


a^^{s,0/d£, 



( 22 ) 


taking ^ = 7TeV and = 8,13TeV. The out¬ 
come is plotted in the upper panel of Fig. To first 
approximation, towards large ^i _2 the effective scaling 
exponent is aeff ~ 0.35 and independent of the lep¬ 
ton charge. In more detail, the scaling exponent is not 
exactly constant but some variation is visible which re¬ 
flects the fact that the PDFs do not follow a pure power 
law, especially when x is not very small (at small ^ 1 , 2 )- 
The scaling exponent for tends to have more slope 
and to be somewhat larger than that of especially 
at small ^ 1^2 which corresponds to midrapidity. This 
can be explained by the slightly steeper slope of the 
u distribution in comparison to d distribution (see the 
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Fig. 3 Lepton charge asymmetry in p-p = 1.96 TeV), p— 
p(^ = 1.96,7,8TeV), p-Pb (^5 = 5.02TeV) and Pb-Pb 
= 2.76 TeV) collisions, for y > 0 (upper panel) and y <0 
(lower panel). 


lower panel of Fig. and the fact that £~ production 
tends to be sensitive to somewhat larger values of x 
in the small-a; side. The latter follows from the factors 
(1 ± a/ 1 — that multiply PDFs in Eq. (j^. 

These, in turn, originate from the parity non-conserving 
W couplings to quarks and leptons. The lower panel in 
Fig. [2] compares the extracted exponent aeff ~ 0.35 to 
the CTIONLO sea-quark PDFs. Evidently, there is a 
good correspondence between the scaling exponent a 
and the behaviour of the small-cc quark PDFs. We can 
conclude that despite the complex higher-order QCD 
calculations, the centre-of-mass dependence of the cross 
sections being discussed can be essentially captured by 
a simple power law. 


Let us now discuss Eq. (21) and whether the na¬ 
ture of the hadronic projectile or nucleus probed at 
small X really disappears as conjectured. To this end we 


have computed the lepton charge asymmetry (again, at 
NLO accuracy) in various collision systems at centre-of- 
mass energies that correspond to existing Tevatron and 
LHC data. The results are shown in Fig. At y 0, 
the curves corresponding to p-p, p-Pb and p-p tend to 
unite, whereas in the opposite direction, y ^ 0, p-Pb 
and Pb-Pb become approximately the same. Thus, as 
far as theoretical NLO expectations are concerned, the 


scaling law of Eq. (21) turns out to be a very good 


approximation, though not perfect. The largest devia¬ 
tions in Fig. are seen in the case of p-p at the Teva¬ 
tron energy, ^/s = 1.96 TeV. There, the probed values 
of X for p are not small enough and especially the as¬ 
sumption of charge-conjugation symmetric quark dis¬ 
tributions, Eq. (0, is not particularly accurate until 
almost the end of phase space (e.g. = 1 corresponds 

to X 2 « M^/s « 0.002). The p-p curve at the same 
center-of-mass energy unites with the rest already at 
lower ^ 1 . 

At small fixed value of the lepton charge asymme¬ 
try in p-p collisions tends to decrease towards increas¬ 
ing centre-of-mass energies. This can be interpreted in 
terms of slightly different scaling exponent for and 
i~ production (see Fig. [^. Denoting the scaling expo¬ 
nent for production by and the difference by 
A = a~ — 




to first approximation, 

Hi,H2/ 


(23) 


1 - 


Z\ 

Since Z\ > 0 


C^’“MAC)]'|log(^) +0{A^). 


we have a condition 




(24) 


which explains the decreasing trend of lepton charge 
asymmetries in p-p collisions towards higher centre-of- 
mass energies at fixed, small 


4 Data and predictions 

4.1 Comparison with existing data 

The currently most accurate experimental measurements 
for inclusive W production from Tevatron and LHC ex¬ 
periments are summarized in Table A direct com¬ 
parison of various measurements is complicated by the 
kinematic cuts for lepton pT, missing transverse energy 
.^X) s-iid transverse mass mx of the neutrino-lepton sys¬ 
tem, which vary among the experiments and have to be 
accounted for. Here, we have chosen to “correct” the 
data to pt > 25 GeV (the default cut in CMS measure¬ 
ments) by MCFM evaluating the observables first with 
the true cuts shown in Table then with px > 25 GeV 
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Table 1 The experimental data sets. 


Experiment 

System 


DO 

p-p 

1.96 TeV 

ATLAS 

Pb-Pb 

2.76 TeV 

GMS 

p-Pb 

5.02 TeV 

ALICE 

p-Pb 

5.02 TeV 

CMS 

p-p 

7 TeV 

ATLAS 

p-p 

7 TeV 

LHCb 

p-p 

7 TeV 

LHCb 

p-p 

8 TeV 

CMS 

p-p 

8 TeV 


kinematic cuts Ref. 


PT > 25 GeV, > 25 GeV [5] 

PT > 25 GeV, > 25 GeV, rm > 40 GeV [22] 

PT > 25 GeV [20] 

PT > 10 GeV ED 

PT > 25 GeV [m 

PT > 20 GeV, > 25 GeV, rm > 40 GeV [49] 

PT > 20 GeV ED 

PT > 20 GeV ESI 

PT > 25 GeV ESI 


-Q 

D- 




n ^ 1 11111 

p‘t > 25 GeV 

. 

j @ i 


y>0 




i.. 

i: 




k, 




k: 




i.. 


□ ATLAS, p-p, 7TeV 

i.. 


♦ LHCb, 

p-p, 7TeV 

Hi ' 


X LHCb, 

p-p, 8TeV 


- 

■ CMS, 

p-p, 8TeV 



• CMS, 

p-Pb, 5.02TeV 

X 


* ALICE, p-Pb, 5.02TeV 

!♦... 


CTIONLO, p-Pb. 5.02TeV 

_^^^_i ■ ■ ■ i 1_ 




10 '^ 10 * 1 


-Q 

& 




-C 

. a 



0.5 

0.45 

0.4 

0.35 

0.3 

0.25 

0.2 

0.15 

0.1 

0.05 

0.0 


f I' 




pf > 25 GeV 
y>0 




'"'i-s 



□ ATLAS, 


: 


p-p, 7TeV ' 



♦ LHCb, 

p-p, 7TeV 



X LHCb, 

p-p, 8TeV 

ix 


■ CMS, 

p-p, 8TeV 



• CMS, 

p-Pb, 5.02TeV 


_ 

* ALICE, 

p-Pb, 5.02TeV 



CTIONLO, p-Pb. 5.02TeV 

■■u 



10* 



P., 



0.5 

0.4 

0.3 

0.2 

0.1 

0.0 

10 ^ 10 ‘ 1 
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p't > 25 GeV 
y <0 


□ ATLAS, Pb-Pb, 2.76TeV 

• CMS, p-Pb, 5.02TeV 

* ALICE, p-Pb, 5.02TeV 
CTIONLO, p-Pb, 5.02TeV 


Fig. 4 Absolute spectra of charged leptons (upper panels for £+, lower panels for £~) in p-p (y^ = 7,8TeV) and p-Pb 
(y's = 5.02TeV) collisions for y > 0 (left-hand panles), and in Pb-Pb (y's = 2.76TeV) and p-Pb (y/s = 5.02TeV) collisions 
for p < 0. The data has been scaled by (s/GeV^)“'^-^°. 


and taking the ratio (absolute cross sections) or differ¬ 
ence (charge asymmetry). We stress that if the kine¬ 
matic cuts were the same in all experiments, this step 
would be unnecessary. The available absolute cross sec¬ 
tions are compared in Fig. The p-p and p-Pb data 
are plotted together at forward rapidity (left-hand pan¬ 
els) and Pb-Pb and p-Pb data together at backward 


rapidity (right-hand panels). In these plots, the data 
has been scaled by a factor (s/GeV^)““, where a con¬ 
stant value a = 0.4 has been used for the scaling expo¬ 
nent as a compromise between the expected exponent 
at small and large see Fig. Keeping in mind the 
“non-constantness” of the scaling exponent and that at 
forward (backward) direction the p-Pb (Pb-Pb) data 

























































7 


(U 


c 

o 



> 25GeV 

r; 2/ > 0 

M 

' • : 

^..V:. 

X DOp-p, 1.96TeV 

i : 

Z : 

*■: 

• CMS p-p, 7TeV 


° CMS p-p, 8TeV 

i 

* LHCb p-p, 7TeV 

- 

♦ LHCb p-p, 8TeV 

b - 

■ CMS p-Pb, 5.02TeV 

li 

* ALICE p-Pb, 5.02TeV 


CTIONLO p-Pb, 5.02TeV 

- 

_ ^^^^_ 1 1 1 1 1_ ^_ 

- 


10 '^ 10 ‘‘ 1 
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(U 

£? 




Pt 

> 25GeV 



y <0 


I S:.., 

..T 



[ T 



- 



- 


CMS p-Pb, 5.02TeV | i 

- 


ALICE p-Pb, 5.02TeV 


- 

ATLAS Pb-Pb, 2.76TeV 

- 

- 

CTIONLO p-Pb, 5.02TeV 

- 


. 


I":, , , 


10 '^ 10 "' 1 
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Fig. 5 Lepton charge asymmetry in p-p(\/i = 1.96 TeV), p-p (\/s = 7, 8TeV), p-Pb (y'i = 5.02 TeV) and Pb-Pb (y's = 
2.76 TeV) collisions. The dotted curve is to guide the eye and corresponds to at = 5.02 TeV. 


are presumably affected by small-cc shadowing in com¬ 
parison to p-p (p-Pb), an exact match with p-p (p-Pb) 
is not expected. Nevertheless, there is clearly a rough 
correspondence between the data from different colli¬ 
sion systems and different y/s. 

The data for lepton charge asymmetries Ci are com¬ 
piled in Fig.[^ We note that some experimental uncer¬ 
tainties, luminosity above all, cancel in the measure¬ 
ment of the lepton charge asymmetries as compared to 
absolute cross sections. As previously, the data from 


p-p, p-p, and Pb-Pb collisions are plotted only in the 
direction where they are supposed to merge with p-Pb 
data. To a very good approximation, the experimental 
data indeed line up to the same underlying curve which 
corresponds to the charge asymmetry in p-Pb collisions. 
Two CMS p-Pb data points at negative rapidities ap¬ 
pear to lie below the NLO predictions and could poten¬ 
tially require additional nuclear modifications in PDFs 
(as also pointed ont in Ref. [10]). However, the ATLAS 
Pb-Pb data shows no sign of such a disagreement with 
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CTIONLO p-Pb, 5.02TeV 


V X DO 

p-p, L96TeV 

k. 

JT • CMS 

p-p, 7TeV 

i 

j * CMS 

p-p, 8TeV 


ll * LHCb 

p-p, 7TeV 


♦ LHCb 

p-p, 8TeV 

1 

■ CMS 

p-Pb, 5.02TeV 


^ ALICE 

p-Pb, 5.02TeV 


n ATLAS 

Pb-Pb, 2.76TeV 



12 3 4 

l/ref(v^ = 5-02TeV) 


Fig. 6 The world data on lepton charge asymmetry as a function of taking y'slTf = 5.02 TeV. 


the theory at those values of rapidity indicating that 
there appears to be some tension between these two 
data sets and that the both data sets cannot be opti¬ 
mally reproduced with the same set of (nuclear) PDFs. 

We can also compress all the data into a single plot. 
This is done by choosing a certain reference centre-of- 
mass energy (we take .^Sref = 5.02 TeV) and plot¬ 
ting the data as a function of variable 

j/ref = y± ^log—, 2/^0, (25) 

2 s 

such that 

G(y,vG) = G(2/ref, 2 /> 0, (26) 

6(2/, Vs) = 6(2/ref, Gs6f), 2/ < 0. 

Such a plot is shown in Fig.|^ In order to keep the plot 
readable Pb-Pb data is plotted only at 2 / < 0, and p-p, 
p-p data is plotted only at y > 0. 


4.2 Cross-section ratios 

In Section [01 we noted that ratios of cross-sections at 
two nearby ^/s at fixed values of scaling variable 6.2 
could become less prone to large-cc PDF uncertainties 
in comparison to taking the ratios at fixed rapidity. To 


investigate this statement quantitatively, we have com¬ 
puted (p-p collisions, NLO precision) ratios 


^N/V/^(2/ref) 




{V^)/dyrei _ / 

da^^ {yfs)/dy^ei \ / 

dV (vV)/(i2/ref _ / 

da^~{^/s)/dy,ef \ Vs J 


(27) 


(28) 


(29) 


where the prediction from scaling laws are also indi¬ 
cated. For comparison we evaluate the same ratios also 
at fixed rapidity (instead of hxed i/ref)- We have used 
PDF4LHC15_30NL0 set of PDFs (SO] available from the 
LHAPDF libarary [^. This is a hybrid set that com¬ 
bines [S2] information from independent PDF hts (CT14 
[53], MMHT14 [54j, NNPDF3.0 [55|) thereby giving a 
better idea of the uncertainties than when sticking to a 
one particular PDF provider. 

The results are shown in Fig. 0 (v^ = 8 TeV, v 6 = 
7TeV), Fig.0(vG'= 13TeV, v 6 = 8 TeV), and Fig.0 
{VV = 14TeV, y/s = 13TeV). The histograms in red 
indicate the outcome when the ratios are taken at fixed 
rapidity intervals and the green ones correspond to mak- 
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Fig. 7 Ratios of (left) and l~ (middle) spectra computed at y's = 8TeV and y's = 7TeV center-of-mass energies. In red 
color are the results binned in lepton rapidity y, and in green the results binned in taking .y/s^ef = 7TeV. The dashed 
lines indicate the prediction of scaling law Eq. The right-hand panel shows the double ratio of Eq. (|29[|. 




!;,!/,..t(ViW = 8TeV) 


!/. = 8TeV) 


Fig. 8 As Fig. but using = 13TeV and y/s = 8TeV. 
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ing the ratios at fixed y^ei (equivalent to fixed ^i). One 
can observe that in the case of W“ production and the 
double ratio the PDF uncertainties indeed tend to can¬ 
cel out better when the ratios are taken at fixed y^ef- 
For production it appears that there is no definite 
advantage (in the sense that PDF uncertainties would 


decrease) in binning as a function of yref- We attribute 
this to the fact that in the case of W+, the integrand (in 
Eq. @) in X is broader for production than what 
it is for W“ production and the PDF uncertainties do 
not cancel as effectively. 
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The LHCb collaboration has recently reported m 
ratios similar to ones discussed here (though integrated 
over the rapidity interval 2 < y < 4.5), and has ob¬ 
served some deviations between the measurements and 
NLO calculations. Our results suggest that by making 
the rapidity intervals equal in y^ef, the PDF uncertain¬ 
ties especially in the double ratio can be suppressed and 
the significance of the measurement thereby increased]^ 

5 Summary 

We have discussed the scaling properties of inclusive 
charged leptons from decays of W bosons created in 
hadronic collisions. Based on the leading-order estimate, 
we have found that the ^/s dependence of cross sections 
in forward/backward directions at fixed value of scal¬ 
ing variable ^ 1^2 = (-^w/-\/s)e^^ should approximately 
obey a one-parameter power law, in which the scal¬ 
ing exponent is approximately independent of the lep¬ 
ton charge and reflects the slope of the small-cc PDFs. 
Consequently, the lepton charge asymmetries at dif¬ 
ferent centre-of-mass energies are predicted to be ap¬ 
proximately same at fixed ^ 1 ^ 2 - Moreover, lepton charge 
asymmetries in different collision systems are related: at 
large positive (negative) y the lepton charge asymmetry 
depends effectively only on the nature of the forward- 
(backward-) going nucleon or nucleus. A comparison 
with the experimental data from LHC and Tevatron 
confirms that the derived scaling laws are indeed able 
to capture very well the behaviour of the data. 

While these scaling laws by no means serve as a re¬ 
placement for accurate (NLO and beyond) calculations, 
the possibility of a direct comparison of various data 
should be useful in e.g. checking the mutual compati¬ 
bility since by fixing (^ 2 ) one forces the PDFs to be 
sampled at approximately the same regions of Xi (X 2 ) 
independently of ^/s. This, as we demonstrated, can in 
turn be taken advantage of by reducing PDF uncertain¬ 
ties in ratios of cross sections measured at different y's. 
This could increase the sensitivity of the experiments 
e.g. to possible contributions from physics beyond the 
Standard Model. 
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